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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

THEORETICAL CALCULATIONS OF THE STABILITY DERIVATIVES 
AT SUPERSONIC SPEEDS FOR A HIGH-SPEED 
AIRPLANE CONFIGURATION 
By Kenneth Margolis and Percy J. Bobbitt 

SUMMARY 


Theoretical calculations of the stability derivatives at supersonic 
speeds for a high-speed airplane configuration are presented. The range 
of Mach numbers considered includes those cases for which the wing and 
tail surfaces have subsonic leading edges, supersonic leading edges, and 
combinations of the two. The methods of analysis are discussed; these 
represent detailed consideration of the important effects and contribu- 
tions of the various airplane components utilizing available theories 
and theoretically justifiable and reasonable approximations. 

Results of the investigation, which are applicable to motions 
involving small angles of attack and sideslip and low rates of rolling, 
pitching, and yawing, are presented in the form of graphs illustrating 
the variation of the stability derivatives with Mach number. 


INTRODUCTION 


The prediction of the stability of complete airplane and missile 
configurations requires a knowledge of the aerodynamic forces and moments 
acting on all component surfaces of the airframe and the rates of change 
of these forces and moments with respect to the various attitudes , veloc- 
ities, and accelerations (i.e., stability derivatives). Design of air- 
planes and missiles for automatic stabilization has accentuated the need 
for this information while the design is still in its initial stages. 

Thus the problem of estimating by theoretical means the forces, moments, 
and resulting stability derivatives of complete airplane configurations 
at supersonic speeds is one of considerable interest and importance. 

This is particularly true at the present time in view of the difficulties 
encountered in determining these quantities experimentally for certain 
motions. While there is considerable information available for isolated 
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aircraft components there is relatively little knowledge, either of a 
theoretical or experimental nature, concerning complete configurations. 

As a result, it is extremely difficult to estimate with accuracy the 
performance and stability characteristics of present-day missiles and 
airplanes. 

During the past few years (and at present), considerable research 
effort has been directed towards the theoretical determination of wing- 
body and body-tail interference effects and the effects of induced flow 
on the loading characteristics of the component surfaces of complete con- 
figurations. Accurate estimates of these mutual and induced interference 
effects as well as additional studies related to isolated components and 
components in combination are required before accurate theoretical esti- 
mates of stability derivatives for arbitrary complete configurations can 
be accomplished. 

The present paper results from an investigation which concerned 
primarily the problems involved in adapting, modifying, and extending 
the available theories as Required in order to estimate reliably the 
stability derivatives at supersonic speeds for arbitrary complete air- 
plane configurations. For purposes of application, a high-speed 
interceptor -type aircraft was chosen. It is primarily this phase of the 
investigation, that is, the estimation of stability derivatives for the 
specific configuration, which is reported herein. The methods used are 
described; particularly detailed emphasis is given the case of steady 
rolling. Results of the analysis are presented in the form of graphs 
illustrating the variation of the stability derivatives with Mach number. 


SYMBOLS 


*,y,z 

Cartesian coordinates 



v,w 

perturbation velocities along 
tively; negative values for 

y 

w 

and z axes, respec- 
indicate downwash 


angular velocities about x, 
respectively 

y. 

and z axes, 

a 

angle of attack 




d rate of change of angle of attack wit 1 ' respect to time 

P angle of sideslip 

V flight velocity 
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B = \ (Mach number)^ - 1 


P 

0 


'P 


density of air 

coefficient of pressure difference between upper and 
lower surfaces , positive in the sense of lift 

velocity potential evaluated on upper surface 


circulation,, ^ 


>T.E. 


L.E. 


AC p dx 


S 

A 

b 

b vt 


wing area including portion masked by body 
aspect ratio of wing 
wing span 
vertical-tail span 


c 

l 

A-vt 

C 

K = -JtH 

BA^ - 2BC 


mean aerodynamic chord of wing 

distance between center of gravity and center of pres- 
sure of sideslipping vertical tail 

aspect ratio of vertical tail 
leading-edge slope of vertical tail 


a distance from roll axis to wing-body juncture 

x^y-L Cartesian coordinates used in analysis for rolling motion 


X = x - xp 

Y = y - yp 

i variable index, used in summations and as subscripts 

Xj_ = x - x-j_ 

Yi = y - yi 
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X^Y^z 

m 

k,k x ,k 2 

L 
Y 
L 1 
N 
M 


Xj_, Yj_, and z, respectively , made nondimens ional with 
respect to wing semispan 

slope of lifting line (absolute value) 

constants 

limits of integration 

lift 

side force 
rolling moment 
yawing moment 
pitching moment 


lift coefficient, 


|pV 2 S 


% 


-'n 


“'m 






side-force coefficient, 


§pV 2 S 


rolling-moment coefficient. 


L 1 


|pV 2 Sb 


yawing -moment coefficient. 


N 


§pV 2 Sb 


pitching-moment coefficient, 


M 


|pV 2 Sc 


dC L j 


,5a 


cl — >0 
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SC l\ 

a ^ 

. 2V/r-.->0 

Subscripts : 

t tail contribution 

t,o tail contribution assuming zero end plate 

All angles are measured in radians unless indicated otherwise. 



SCOPE 


The geometry of the configuration under study is shown in figure 1 
and the detailed characteristics of the component surfaces are given in 
table I. Considered in the analysis are the following motions: constant 

angle of attack, steady rolling, steady pitching, constant vertical accel- 
eration, constant sideslip, and steady yawing. The following stability 
derivatives are obtained: C-^, C^, C Yp , C np , C Ip , C,^, C,^, C Yp , 

c np> c Y r ’ c n r > and C V 

All airframe components are considered rigid; wing and tail surfaces 
are assumed to have zero thickness. The wing, horizontal tail, and body 
are at zero geometric angle of attack for all motions other than angle of 
attack. The analysis is carried out within the framework of the linearized 
theory, the results being valid for small angles of attack and sideslip and 
for low rates of rolling, pitching, and yawing. In addition, modifica- 
tions and assumptions are introduced into the analysis as needed in order 
to cope with the problems encountered in estimating the derivatives for 
this particular configuration. Although the final results are by no 
means based completely on theoretically rigorous procedures, they do 
represent detailed consideration of all inqportant effects and contribu- 
tions utilizing theoretically justifiable concepts. 

Calculations are carried out for four Mach numbers: {2, 2, 2.5, and 

It may be noted at this point that the range of Mach numbers considered 
includes the cases where the wing and tail surfaces have subsonic leading 
edges, supersonic leading edges, and combinations of the two. 

Results for the stability derivatives are .presented relative to a 
body system of axes maintaining, however, the usual stability convention 
for denoting positive forces and moments. (See fig. 2.) Radian measure 
is used throughout for the derivatives . 
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STABILITY DERIVATIVES DUE TO ANGLE OF ATTACK 
General Considerations 


The body, wing, and horizontal tail are all assumed to be at a small 
angle of attack. The lift on the isolated body is concentrated mainly 
over the forward portion; the exposed panels of both the wing and hori- 
zontal tail in addition to carrying their own lift are affected by the 
presence of the body and in turn induce some lift on the body. There is 
also, of course, the downwash effect from the wing which decreases the 
lift effectiveness of the horizontal tail. Fortunately, for the angle- 
of-attack case there is, relatively speaking, considerable published 
information to be drawn upon for purposes of estimating the various con- 
tributions. The isolated wing and tail surfaces, both of which are modi- 
fied delta plan forms, are treated in references 1, 2, and 3> wing-body 
and body -wing interference effects may be estimated most conveniently and 
apparently quite satisfactorily for circular bodies of revolution by use 
of references *1-, and 6. The effect of the isolated body is difficult 
to determine accurately because of the cross-sectional deviations of the 
airplane fuselage from the theoretically treated bodies of revolution. 
However, slender -body theory for bodies of revolution in conjunction with 
some results obtained in references 7 and 8 for other types of bodies 
enable a qualitative evaluation of the effects of the isolated body and 
its influence on the mutual interference effects between body and wing or 
tail surfaces. The induced flow from the wing has been treated in some 
detail by various investigators (see, for example, refs. 9 to 12) and 
hence the net contribution of the horizontal tail may be readily calcu- 
lated. It appears, therefore, that by judicious use of the available 
information theoretically justifiable and reliable estimates of the 
derivatives Ci^ and may be anticipated. 


Detailed Considerations and Results 

The lift-curve slope and the pitching moment due to angle of 

attack ^i.e., stability derivative for the isolated wing and iso- 

lated horizontal tail are directly obtainable from reference 1 for the 
subsonic ieading-edge cases and from references 2 and 3 for the supersonic 
leading-edge cases. Although a number of methods were investigated in order 
to estimate the interference effects between wing and fuselage (and between 
horizontal tail and fuselage), the results obtained agreed in general with 
those predicted by use of references k to 6, hence the details of the 
calculations will not be presented. The lift on the isolated fuselage 
was calculated to be slightly greater than that predicted by slender - 
body theory. Inasmuch as slender-body theory fails to predict any lift 
on the cylindrical portion of a cone-cylinder fuselage combination, it 
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is felt that an estimate based in part on the results given in refer- 
ence 7 is more realistic. The downwash in the vicinity of the horizon- 
tal tail was calculated using the procedures indicated in references 9 
and 10. The assumption (which for this particular configuration is quite 
reasonable) is made that the vortex sheet remains flat in the vicinity of 
the tail. Values of the downwash were obtained at several stations behind 
the apex of the horizontal tail and then averaged over the horizontal tail 
surface to form an "effective" downwash velocity. Figure 3 presents the 
variation of the downwash factor w/aV with Mach number. The corresponding 
induced angle of attack at the tail was then added to the geometric angle 
of attack of the tail and the resulting lift (neglecting the presence of 
the body) was calculated using references 1 and 2. Interference effects 
between horizontal tail and fuselage were calculated in a manner similar 
to that used for the wing-fuselage interference effects. The net result 
of adding the lift contributions of all components taking into account 
the calculated interference and induced flow effects, that is, the estimated 
lift-curve slope of the complete airplane configuration, is presented in 
figure 4 - together with some breakdown results. The ordinate is the lift- 
curve slope based on a reference area S that includes both the 

area of the exposed wing panels and the masked (by the fuselage) portion 
of the wing extended to the fuselage center line. As mentioned previ- 
ously (see section entitled "SC0PE n ), actual calculations were carried 
out for four values of the abscissa Mach number (\[2~, 2, 2.5; and 3) and 
the final curve faired accordingly. 

Calculation of the pitching moment due to angle of attack is a bit 
more difficult inasmuch as detailed information for the incremental center - 
of -pres sure shifts or incremental pitching moments due to the various com- 
ponent effects discussed previously ^for C is lacking. Reference 3 

considers the interference effects between wing and fuselage (and between 
horizontal tail and fuselage). Calculations of the other effects were 
made by first assuming center-of -pressure locations (based on approxi- 
mate lift distributions due to each effect or component) and then evalu- 
ating the significance of varying the assumed center-of-pressure locations 
through a reasonable range. Inasmuch as the pitching moment is signifi- 
cantly different from zero, this method of "averaging" the results should 
yield satisfactory estimates. The resulting pitching -moment derivative 
for the complete configuration expressed in derivative form is presented 
in figure 5; some breakdown results are also shown. The stability deriv- 
ative Cj^ is based on the reference area S and the mean aerodynamic 

chord c. Moments are taken about a center-of-gravity location assumed 
to be at 32 percent of the mean aerodynamic chord. 
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STABILITY DERIVATIVES DUE TO PITCHING 
General Considerations 


Aside from calculations for isolated wings (e.g., refs. 1, 3 > 
and l4) there is relatively little theoretical information available con- 
cerning the pitching-moment derivatives C m ^ and C m ^. Some results 

for various types of isolated bodies of revolution are given in refer- 
ences 15 and l 6 . Reference 17 considers time -dependent downwash at the 
tail, which requires knowledge of the downwash due to steady pitching as 
well as that due to constant angle of attack. Calculations for angle -of - 
attack downwash are available from a number of sources (see section 
regarding angle of attack); data on steady -pitching downwash are lacking 
although rather simplified approximations have been used to advantage in 
the past. More detailed calculations of the downwash due to steady 
pitching and the downwash due to linear angle -of -attack variation with 
time (i.e., constant a) are in progress. Although subject to restric- 
tions on the permissible configurations, reference 18 is quite useful in 
determining wing-body effects. Thus, the estimation of the derivatives C^ 

and C^ for an arbitrary complete configuration is dependent on essen- 
tially the material and references given above. Although the present 
wealth of published or near-published theoretical studies for pitching 
motion is in no way comparable to that for angle of attack, it has been 
shown in an unpublished analysis to be adequate for the prediction of the 
level or magnitude of Cm^ + C^ for a number of missile configurations. 

In view of this evidence there is reason to believe that the theoretical 
estimates of the pitching derivatives for the present configuration are 
fairly realistic in spite of the fact that its geometric variations from 
theoretically treated cases are significant in some respects. 


Detailed Considerations and Results 
The derivatives and for the isolated wing and horizontal 

tail are directly obtainable from reference 1 for the subsonic-leading- 
edge cases. For the case of a supersonic leading edge C m ^ is given in 

reference 3 and C^ in an unpublished analysis. ^Actually the wing and 

horizontal tail are so near to being delta plan forms that the formula 
given in reference 14 for may be used without significant error. ^ 

In calculating the wing -body effects by use of reference 18, the wing 
trailing edge was modified slightly to form a delta wing with the same 
area as the actual wing. The tail contributions to the derivatives C m ^ 
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and were calculated from simplified approximations based in part 

on reference 17. The tail contribution to C™ was obtained by adding 


to the isolated-tail value a correction for the wing downwash based on 
equation (21 ) of reference 17* Inasmuch as the horizontal tail of the 
airplane configuration is only a short distance behind the wing, the 
degree of accuracy of this approximation required consideration. For 

a subsonic -edge delta wing at Mach number \f 2, values of w^ based on 


more exact procedures were compared with those used herein 
correction to 


The resulting- 


Cjh based on the more exact approach did not result in a 
significant change. The tail contribution to was calculated from 


an unpublished analysis which was based on a simplified approximation to 
the method presented in reference 17. The downwash due to & also 

checked fairly well at M = ^2 with more exact calculations. Summation 
of the results obtained for the various isolated components, components 
in combination, and induced-flow effects yielded the curve shown in fig- 
ure 6 for the damping in pitch of the complete configuration. 

A partial breakdown of the results is also indicated. It is of interest 
to point out that the contribution from the & term alone is unstable 
for this particular configuration but fortunately is small in comparison 
with the damping effect due to the steady pitching component. 


STABILITY DERIVATIVES DUE TO STEADY ROLLING 
General Considerations 


In order to estimate with a reasonable degree of accuracy the sta- 
bility derivatives associated with steady rolling motion Cy^, 

and C n ^ for a complete configuration, a prime requirement is to have 


reliable knowledge of induced-flow effects. Unfortunately, there is a 
lack of information with regard to this aspect. Results for isolated 
wings (and horizontal tails) of various plan forms are available in a 
number of papers (e.g., refs. 19 to 23); wing-body and tail-body inter- 
ference effects for restricted types of configurations are treated in 
references 2k to 2 6. Various types of tail configurations are considered 
in references 27 to 29. An attempt to take induced-flow phenomena into 
account in estimating the derivative Cj for a missxle configuration 


was made in reference JO, but the results were not too general in nature 
nor was the accuracy of the procedure involved satisfactory. Thus, 
although for the previous motions treated (angle of attack and pitching) 
experience had shown the methods employed therein in estimating the 
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derivatives to be reliable, such is not the situation for the case of 
rolling. Consequently , detailed attention will be given to the deriva- 
tion of expressions relating to induced-flow calculations as well as to 
the estimation of the contributions due to the pertinent airplane com- 
ponents and components in combination. 


Detailed Considerations and Results 

Wing loading . - In order to obtain' reliable force and moment deriva- 
tives for a wing -body- tail configuration in roll, the effect of the flow 
field on the tail must be considered. Since the flow field is determined 
by the span loading, it is desirous that it (the span loading) be calcu- 
lated as accurately as possible. Large changes in the spanwise loading 
near the root of a wing which might be occasioned by the presence of a 
body or end plate may cause a noticeable difference in the flow field, 
yet affect the damping of the wing itself only slightly. For this 
reason it was considered necessary to modify the isolated-wing loading 
to take into account the presence of the body . 

In addition to this ’’end-plate" effect, the presence of the body 
causes a change in the position of the leading Mach lines on the wing 
(as compared to the isolated-wing case). Whereas for the isolated-wing 
situation the leading Mach lines emanate from the wing apex, the leading 
Mach lines for the wing-body combination originate at the intersections 
of the wing leading edge with the body. 

In order to approximate the effect on the wing loading of super- 
posing a body (and the corresponding change in Mach line position), a 
rather straight-forward procedure has been followed. First, the span 
loading in roll over a wing composed of the two wing panels exterior to 
the body and joined together at their roots was determined by use of 
references 31 and 32. Then the two panels were separated and placed on 
the body. Since the previously determined span loading for the rolling 
wing was predicated on the wing rolling about its axis of symmetry, that 
is, the half -wings rolling about their root chords, an additional loading 
was added to compensate for the greater distance of each point on the 
wing from the roll axis. (The angle -of -attack distribution which gives 
rise to this additional loading is constant and equal to the quantity pa/V.) 
This procedure did not predict the loading on the body. In order to com- 
plete this portion of the loading curve, a straight line was drawn between 
the finite values at the wing -body juncture and zero at the roll axis. 

The resultant span loadings for the wing-body configuration at Mach num- 
bers of {2, 2, 2.3, and 3 are shown in figure 7. For purposes of com- 
parison, figure 8 presents the loading at Mach number ^2 obtained by the 
present method and the corresponding loading for a fictitious isolated 
wing formed by extending the leading and trailing edges of the exposed 
wing panels to the center line. It will be noted that the loading at the 
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root of the wing predicted by the present method is slightly higher than 
that obtained for the isolated wing. As the distance from the root along 
the span increases , the difference between these two loadings decreases. 
Further calculations and comparisons at the higher Mach numbers indicated 
insignificant differences between the modified loading (present method) 
and isolated wing loading at all spanwise stations. 

As expected, it was found that the improved or modified loadings 
yielded wing-body force and moment derivatives that differed only slightly 
from those obtained by use of isolated-wing loadings. However, it appeared 
that the modified loading curves would produce a noticeable change in the 
flow-field calculations, especially at the lower Mach numbers, and hence 
they have been utilized in all such calculations. 

Flow-field theory . - There are many papers which present general for- 
mulas for the potential in space due to a thin lifting surface. In most 
cases these formulas are extremely difficult to evaluate. It becomes 
necessary, therefore, to utilize some approximation which may be more 
easily evaluated and which lends itself to numerical calculations. 

Reference 9 indicates that a lifting line can be used as a very good 
approximation for most downwash problems. A similar statement regarding 
sidewash problems may not be made since there are no exact expressions 
available for the sidewash behind rolling wings. It might be noted, how- 
ever, that an approach somewhat like the above was used in reference 33 
for calculating sidewash behind rolling wings traveling at subsonic speeds 
and very good agreement with experimental results was realized. 

Taken together, references ^k , 9, and 12 represent a fairly thorough 
study of the lifting-line approximation, especially with regard to down- 
wash problems, and show that bent lifting lines will probably give the 
best results for swept and triangular wings. The shape of the bent 
lifting line most often used in these papers and also used herein con- 
sists of a pair of straight lines connecting the midpoints of the root 
and tip chords of the wing. The potential due to a slanted lifting line 
of this type may be obtained from reference ^k (certain errors present 
therein having been corrected in the following equation): 



JL_ 

2 * 



dr ( y l) tan -l z\X 2 - B 2 (y 2 + z 2 ) 


d yi 


z^ Y 

YX - — - ±_ 


dyi 


( 1 ) 
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where the equation of the lifting line is 


X 1 


71 + k 
m 


(See fig. 9 for pertinent symbols, etc.) 
tion (l) becomes 


When 


ar (yi) 

dyi 


is zero, equa- 


T , -l z 
— tan - 

2 * 


;{x2 - B 2( Y 2 + z 2) 


h 2 


m m 


bn 


( 2 ) 


which might be considered as the potential in space (at a point (x,y,z)) 
of a finite bent vortex of constant strength. A number of these super- 
posed on one another can be used to approximate the potential in space 
of a lifting line with any prescribed lift distribution, for example, 
such as those previously calculated (fig. 7) for the rolling wing -body 
of the configuration (see ref. 3^)* 

From equation (2) the sidewash and downwash due to a bent vortex 
are readily obtained by taking the derivative with respect to y and z 
The following results are obtained: 


r 

£(2X2 • 

m \ 

- B 2 ^ 2 - B 2 z 2 ) - zx(x 2 - B 2 z 2 ) 

2* 

{x 2 - B 2 (y 2 + z 2 ) 

- 4 - tf + z2 ( x2 - B2y2 - B 2 * 2 )] 


(3) 
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Since the loading on a rolling wing is antisymmetrical, the induced 
sidewash velocity from each panel is in the same direction and equal in 
the y - 0 plane (the plane of the vertical tail). For this reason, it 
is necessary to calculate only the sidewash from one panel ('using eq. ( 3 )) 
and double it. To calculate the total downwash at any arbitrary point 
on the horizontal tail requires that the contribution from both panels 
be determined. Hence it is necessary to obtain an expression for 0 Z 

when bent lines have negative slopes. Inasmuch as m and k are 
defined as absolute values, this condition is represented by the fol- 
lowing. equation for the lifting line: 


-yq + k 
X 1 = m 


Appropriate modification of the potential (eq. (2)) and differentiation 
with respect to z yields 

*2 


h l 

( 5 ) 


r 

yP . 

'£(z 2 - 
m \ 

Y 2 ) - B 2 YX 

f 2 

r 

+ 2z 2 ) 

b2y 2 

m 



2 k * 

fx 2 - B 2 

'y 2 + z 2 ' 



z 2 

Y 2 ^ 
m j 

2 

+ z 2 

-(x 2 - b2y 2 - B 2 2 

•1 


Flow-field calculations . - Equation (3) can be utilized to formulate 
an approximate expression for the sidewash due to a series of constant- 
strength, horseshoe vortices spaced along a line so as to represent as 
closely as possible some prescribed span loading. This expression is 


V 


.Y_ r ( y i-n) - r (yi-i) 

60 Atc 


% \2Xi 2 - *1? - B2 Z 2) 

| - zX i (x i 2 _ b2 Z 2] 



( Z 2 Y 

PA ' m 1 

■2\ 2 

r) + * 21 

[xi 2 - B 2 Y ± 2 - B 2 z 2 ) 
(6) 

^2 _ B^^Yi 2 + z2^ 


where Y ± = y - yq and X.j_ = x - x i . The subscript i takes on all 

integral values from 0 to n. Since span load distributions are usually 
determined from prescribed angle-of -attack distributions and since the 
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primary purpose of calculating the flow field behind the wing is to 
determine the induced loads on the tail, it is advantageous to write 
equation (6) in the following form: 


n 

r( y i+i) r ( y i-i) 

fliAx 2 -B 2 Y , 2 -B 2 z 2 ^ 

- zXi^Xi 2 - B 2 z 2 ) 


v/V _ \ 

J?(b/2) 2 p(b/2)i 

mV 1 / 

Hi i=0 

7- _ 52 Y 2\ 2 

^i X i -TT-^-j +z2 

(x-j 2 - B^i 2 - B 2 z 2 ) 
(7) 


From this equation the induced angle of sidewash v/V per unit pb/2V 
is arrived at directly. The bar simply denotes that the quantity is non- 
dimensional with respect to the wing semispan. Following this same pro- 
cedure, the induced angle of downwash per unit pb/2V at any point in 
space behind a rolling wing may be obtained from equations (4) and ( 5 ) as 


0 

r(y i+ i) r(y 1 . 1 ) 

x. 2 

Y.X-5 - - 1 r 

W /V _ _ \ 

p(b/2) 2 p(b/2) 2 

YlXi m V 


iiXi 3 - ^|-(s 2 - Yi 2 ) - B2YiXi(?l 2 + 2Z 2 ) - th _(z 2 + Y ± S ) 


Pb 

2V 


i= n l 


MX-j 2 - b 2 (y., 2 + Z 2 


,j + Z 2 . *i 2 

+ TT + nr, 


+ z^X-j 2 - B 2 Y ± 2 - B 2 z 2 



r (yi+i) r (yi-i) 

— 2 2 - 2 

Yi v i^ 1 ^ f " 2 Yi 2 ^ E 2 Yl v i ^Yi 2 1 °" 2 ^ I ^ ^ (" 2 

+ Yi 2 ) 

^(b/2) 2 p(b/2) 2 

1 i A I t ii 1 d y + in \ z 


(% - C - ¥) + 52 ( J i 2 - B^l 2 - l**) 


i^Xi 2 - b 2 ( 

Xi 2 + 


(8) 


where nq is associated with negative values of y^ and n 2 is asso- 
ciated with positive values of yq. It should be noted that in equa- 
tion (8) the parameters i, nq, and ng are always positive. The 
expression for Xq corresponding to positive values of yq is 


X-r 


yj + k 
m 
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and that corresponding to negative values of y^ is 


x i 


-yi + k 

m 


Seventeen horseshoe vortices on each half -wing were used to approximate 
the span load distribution (expressed in terms of r in the equations). 
(See fig. 10.) Calculations of the induced angles of sidewash and down- 
wash (eqs. ( 7 ) and ( 8 ), respectively) were carried out for a number of 
longitudinal locations along the tail surfaces. For each longitudinal 
location considered, the induced angles of sidewash and downwash were 
obtained at the root, tip, and intermediate spanwise positions of the 
appropriate tail surfaces. By means of an averaging process, resultant 
curves for the spanwise variation of the "effective" induced angles were 
obtained. The results at various Mach numbers for the spanwise variation 
in downwash angle for the horizontal tail are presented in figure 11 ) 
corresponding results for the sidewash angle for the vertical tail are 
presented in figure 12. It should be noted that the chordwise variations 
of the induced angles on the tails were small and hence the curves obtained 
by the averaging process represent fairly accurately the spanwise varia- 
tion at any longitudinal station. It should also be noted that the effect 
of the body on the flow field behind the wing has been assumed to be small 
and therefore neglected. 

Tail loads, forces, and moments .- To determine the induced pressures 
and the corresponding induced loadings on the tail surfaces by the exact 
methods requires detailed knowledge of the induced angles of sidewash and 
downwash at each point of the pertinent tail surface. Furthermore, deri- 
vation of the induced pressures from the induced-angle distributions 
entails considerable difficulties except for the most simple cases. To 
utilize a theoretically rigorous approach in the present investigation 
would obviously necessitate an almost prohibitive amount of work. An 
alternate method of attack on the problem is to determine the desired 
induced loadings on the tail surfaces by utilizing loadings that have 
been previously obtained in references 31 and 52 for isolated wings under- 
going various motions. In order to do this, however, it is necessary that 
the induced-angle distributions for the tail be defined by constant or 
linear (zero value at origin) expressions as is the case, for example, of 
an isolated wing at a constant angle of attack, in steady rolling, or in 
steady pitching. Inasmuch as superposition of these solutions is permis- 
sible, the criterion for enabling the utilization and application of 
isolated-wing loadings to the calculation of induced tail loadings is 
that the induced angles of sidewash and downwash on the tail surfaces be 
approximated by arbitrary straight lines. While it may not be possible 
to represent accurately with straight lines the induced angles of sidewash 
and downwash for all configurations, it does seem that for the majority 
of cases this may be done without incurring a large percentage error in 
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the induced loading. When there is little or no longitudinal variation 
in the induced angles of downwash and side-wash (as has been stated in the 
previous section to be the situation for the configuration under study) 
the induced loading may be built up solely of two load distributions, one 
corresponding to that obtained for a wing at a constant angle of attack 
and the other corresponding to that obtained for a wing in steady rolling. 
This point will be more fully explained in the discussion which follows. 

The straight-line approximation of the average indue ed-downwash and 
induced- sidewash curves for the tail surfaces of the configuration is 
indicated in figures 11 and 12. For convenience, a typical curve has 
been selected (see fig. 13 ) to indicate how the induced-angle distribu- 
tion may be broken up into two components, one component being a hori- 
zontal line defined by 



( 9 ) 


and the other component being a line of nonzero slope passing through 
the origin: 



( 10 ) 


The difference k n - ko is the actual induced-angle distribution 

^ b/2 

(as approximated by a straight line). It is readily seen that the com- 
ponent angle distributions defined by equations (9) and (10) are simi- 
lar to the distributions obtained for isolated wings at a constant angle 
of attack and in steady rolling, respectively. Taking into account the 
constants used to nondimens ionalize the quantities given by equations ( 9 ) 
and (10), the induced load distribution for each component may easily be 
obtained from charts presented in references 31 and 32 . 

The induced load distributions thus far obtained have utilized wing 
loadings that are applicable to one panel of a symmetrical wing. It is 
of course necessary to modify these results so as to account for the 
mutual interference effects between the horizontal and vertical tail sur- 
faces as well as to simulate in some manner the proper end-plate effects 
caused by the presence of the body. Inasmuch as the isolated tail sur- 
faces performing a steady roll have load distributions which must also 
be modified to take into account interference phenomena, it is more 
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convenient to add together the uncorrected induced load distributions to 
the isolated tail load distributions and then correct all load distribu- 
tions simultaneously. The load distribution due to the isolated tail 
surfaces in roll has been obtained directly from references 31 and 32, 
once again considering each tail panel as one-half of a symmetrical wing. 

Inasmuch as there is little information at supersonic speeds con- 
cerning horizontal-vertical tail combinations in roll with the exception 
of "slender" arrangements (ref. 28), the interference effects must be 
evaluated in an approximate manner. Since the tail surfaces of the 
present configuration are highly swept and of moderate aspect ratio, it 
has been found advantageous to modify the previously obtained load distri- 
butions by interference correction factors based on results obtained for 
slender tail arrangements in steady roll (ref. 28) and slender tail 
arrangements at constant sideslip (ref. 35)* 

Before discussing the interference correction factors in detail it 
should be remembered that each induced load distribution was subdivided 
into two component load distributions, that is, an angle-of -attack distri- 
bution (which is directly comparable to a load distribution due to constant 
sideslip for a vertical tail) and a steady-roll distribution. Inasmuch 
as the mutual interference effects are different for the rolling and 
angle-of -attack (or sideslip) cases, each of these load distributions 
must be modified in a different manner. The rolling component of the 
induced load distributions may of course be combined with the isolated- 
rolling-tail load distributions for purposes of applying the interference 
correction factor. 

For the rolling load distributions, interference correction factors 
based on slender-body theory have been used. To obtain these slender- 
body correction factors, the tail system of the configuration under study 
has been approximated by an equivalent tail system composed of half-delta 
panels intersecting in a common chord at the roll axis. The areas of 
both the horizontal and vertical tails in the equivalent system were 
essentially unchanged from those in the original configuration. Actually 
the vertical tail which is by far the more important contributor to the 
rolling moment and the sole contributor to the yawing moment was changed 
only slightly from its original shape. With the ratio of the vertical- 
tail span to the horizontal-tail span known for this equivalent tail 
system, correction factors have been obtained from references 28 and 35 
and applied to the uncorrected load distributions in the following manner: 
Let Lq 1 be the uncorrected rolling load component and L^_ ' be the 

uncorrected angle -of -attack (i.e., sideslip) load component; 
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For the vertical tail. 



L'l of vertical panel of slender inverted T-tail in roll\ 
N 1 


of one panel of a slender triangular wing in roll 


and 



pT 1 

L i 

\ 


L x 'z 

> dz 


X 
1 — 1 

i 

V 

j 

/ 



(XI) 


\ 

of vertical panel of slender inverted T-tail in sideslip 


of one panel of a slender triangular wing at an angle 

of attack 


/ 

( 12 ) 


For the horizontal tail, 

(h' of horizontal panel of slender inverted T-tail in roll\ 


Lo'y iy 


lL t of one panel of a slender triangular wing in roll 


( 13 ) 


The rolling moment contributed by the angle-of -attack component of 
the load distribution on the horizontal tail is proportional to 


J Lx'y d y 


(iM 


There is also an induced loading on the horizontal tail due to the induced 
angle-of -attack component of the load acting on the vertical tail. This 
effect would be in addition to the result given by expression (lif) and has 
been calculated by use of reference 33* 
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The expressions (ll) to (lA) must of course be nondimens ionali zed 
and combined properly to yield tail forces and moments. The tail con- 
tributions to the stability derivatives Cy^ an ^ ^np W then 

readily be obtained by use of the conventional definitions. 


In determining the total forces and moments contributed by the tail 
surfaces, however, only those portions of the loads exterior to the body 

^body boundary defined by — ^ = 0.106 for the case of the vertical taily 

were integrated (see expressions (ll) to (l^i-)). It was assumed that the 
resultant of the forces and moments acting on the body, that is, those 
due to the isolated body in roll, the induced effects, and so forth, was 
very small and hence neglected. 

It is possible to compare the accuracy of some of the approximations 
used with calculations based on rigorous applications of linearized theory. 
Inasmuch as these checks deal with the angle-of-attack component of the 
induced load on the vertical tail (i.e., sideslip effect), the compari- 
sons will be made in a subsequent section of the report dealing with the 
sideslipping motion. 


Specific derivatives .- The value of Cj due to the wing-body com- 
bination is readily obtained from the loadings which were discussed in a 
previous section and presented in figure 7* The addition of the wing- 
body contribution to the total tail contribution yields the damping- in- 
roll derivative C for the complete configuration (fig. 14). It is 

interesting to note that the net tail contribution to the derivative 

was found to be insignificant through the Mach number range considered. 
For this particular configuration, the flow field and other interference 
effects almost entirely nullified the damping of the isolated tail unit. 
Thus the curve shown for the complete configuration is very close to the 
result obtained for the wing-body contribution. 


Inasmuch as the configuration was assumed to be at zero geometric 
angle of attack, there were no suction forces on the wing leading edge 
at Mach numbers for which the edge was subsonic. (At Mach numbers for 
which the wing leading edge is supersonic, there would be no suction 
forces present at any angle of attack. ) Thus the derivatives Cy^ and 

Cnp are com P°sed solely of the tail contributions. The variation of 


the stability derivative Cy^ 


with Mach number for the complete con- 


figuration is shown in figure 15 ; some pertinent breakdown results are 
also included. 


In order to obtain the stability derivative Cnp, each of the previ- 
ously discussed load components acting on the vertical tail was isolated 
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and the center-of -pressure location for each component was determined at 
the four Mach numbers considered. (it is interesting to note that the 
center-of -pressure locations for all components and at all Mach numbers 
did not vary too greatly from the ^A-^oot-chord position.) Each load 
component was then multiplied by the appropriate arm length (measured 
from the center of gravity to the calculated center-of-pressure location) 
and the results added to obtain the total yawing moment. The variation 
of the corresponding derivative C^ with Mach number is presented in 

figure l6. The sidewash effect on the derivative C n ^ ^and 

quite marked; note the reversal in sign near Mach number 1.6. 


%) 


is 


STABILITY DERIVATIVES DUE TO CONSTANT SIDESLIP 
General Considerations 


The information available pertinent to the theoretical calculation 
at supersonic speeds of sideslip effects for complete configurations is 
of relatively scant nature. Results for isolated wings without geometric 
dihedral are available for combined sideslip and angle -of -attack motion 
(refs. 36, 37> and 38) • Wings with geometric dihedral are considered in 
references 39 and 40. For bodies of revolution, results available for 
angle of attack (e.g., refs. 7 and ^l) may of course be applied to the 
sideslip motion. For fuselage shapes with noncircular cross sections 
and abrupt contour changes there is very little information. Some 
results of investigations for various tail configurations in sideslip 
are given in references 33 j 4l, and k2; other information is as yet 
unpublished. Flow- field and component interference effects have not 
received any significant attention to date in the literature. It is 
seen, therefore, that considerable additional research is required before 
accurate estimation of the stability derivatives due to sideslip for 
arbitrary complete configurations can be obtained. In analyzing the 
present configuration the most important effects have been accounted for 
by utilizing the available theory in conjunction with several semiempirical 
approximations. The reasonableness of the approximations is demonstrated 
wherever possible. 


Detailed Considerations and Results 

The forces and moments produced by the vertical tail at an angle of 
sideslip may be obtained by a procedure analogous to that used in treating 
the angle-of -attack component of the induced load on the vertical tail due 
to steady rolling. I^See section for rolling motion entitled ’Tail Loads, 
Forces, and Moments. ) From this angle -of- attack component, forces and 
moments on the sideslipping tail may be obtained in a fairly simple manner. 
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While before the angle of attack was given in terms of pb/2V (see 
eq. (9)) and only a fractional part of the loading due to a fictitious 
angle p was used, now the actual angle of sideslip is (3 and the 
entire loading must of course be considered. It should be noted that 
end-plate effects on the load distribution due to sideslip have also 
been taken into account in treating the angle-of-attack component of 
the induced load due to rolling motion and hence no further considera- 
tion of them is necessary. 


From the loadings just discussed,, contributions of the tail unit to 
the stability derivatives Cy^, C n ^, ma ^ rea dily obtained. 

The accuracy of the end-plate corrections made in obtaining the vertical- 
tail contributions to the sideslip derivatives may be examined in the 
light of some rigorous results available for isolated vertical tails 
with zero and complete end plates. Values for the tail contributions 
obtained herein using wing loadings corrected for end-plate effects 
should, of course, fall between these two limiting cases. 


For isolated vertical tails with supersonic leading edges, both the 
zero- and complete- end-plate solutions for ^Cy^^, ^ , and (plp)t 

are presented in reference 42. For the subsonic-leading-edge condition, 
complete- and zero-end-plate boundaries for ^Cy^^. and may be 

obtained by use of the and Cj^ expressions for isolated wings 

presented in reference 1. (The angle-of-attack solution for a wing is 
comparable to the complete -end-plate case for a vertical tail in side- 
slip.) The complete-end-plate solution for (pl$)t has been derived 

in an unpublished analysis. The zero -end-plate boundaries for the 
derivatives ^Cy ^C n and that are valid at Mach num- 

bers for which the tail leading edge is subsonic may be obtained from 
expressions presented in reference 43* 

The calculated tail contribution (which includes a finite end-plate 
effect) to the stability derivative Cy^ is compared with the zero- and 

complete- end-plate boundaries in figure 17; the reasonableness of the 
approximations used in the calculation is readily apparent. Similar 
comparisons for the tail contributions to the stability derivatives C n p 

and C ^ yielded equally satisfactory results. (it should be mentioned 

that the zero-end-plate boundary for Mach numbers corresponding to sub- 
sonic leading edges was calculated from ref. 43, with the vertical tail 
being slightly modified to form a half-delta of equivalent area.) 
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In calculating the isolated-body effects, it was realized that the 
airplane fuselage would probably have a relatively large side force 
associated with the sideslip motion. (As can be seen from fig. 1, the 
cross sections of the rearward part of the body are almost rectangular 
in shape.) Unfortunately most of the existing information concerning 
bodies is restricted to slender bodies of revolution. No formulas are 
available at present which will predict the forces and moments on yawed 
bodies with noncircular cross sections and abrupt contour changes simi- 
lar to the body considered herein. It becomes necessary therefore to 
approximate the body by an equivalent body of revolution. Obviously 
there are any number of circular body shapes which could conceivably 
represent an "equivalent" body. Thus any theoretical approximation 
based on conventional bodies of revolution would be subject to considera- 
ble uncertainties. In order to minimize the error involved, two estimates 
of the body contribution were calculated; one that was felt to be an 
underestimate and one an overestimate, both being calculated from slender - 
body theory (for example, see ref. 4l). The resulting two values for Cy 


were averaged and this average value was used for the body effect. A 
similar sort of averaging process was used in estimating the longitudinal 
center-of-pressure location in order to obtain the body contribution to 
C n ^. The body contribution to C^p is very small and hence can be neg- 
lected. In connection with these isolated -body effects, it should be 
borne in mind that any inaccuracies present in the body estimates would 
of course affect the total 0 


Y c 


estimate more than the total C 


f-p n P 

total C^p estimates because of the relatively short moment arm involved 

for the body as compared to that of the vertical tail. Variations of the 
derivatives G^, C n p, and C^ with Mach number for the complete con- 
figuration are shown in figures 18, 19 , and 20, respectively; some break- 
down results are also included. 


or 


STABILITY DERIVATIVES DUE TO STEADY YAWING 
General Considerations 


With regard to steady yawing motion, there is an almost complete 
lack of the information required for estimating stability derivatives 
for complete configurations. In fact, this statement is almost equally 
applicable to isolated components. Isolated wings have been treated in 
references 1, 44, 45, and 46 but some of the results obtained therein are 
based on fairly crude approximations and must be considered of only quali- 
tative value. For isolated bodies of revolution the relatively few calcu- 
lations available for steady pitching motion (e.g., refs. 15 and l6) may 
of course be applied to the case of steady yawing. For fuselage shapes 


CONFIDENTIAL 



2k 


CONFIDENTIAL 


NACA RM L53G17 


other than bodies of revolution there is essentially no information. 
Wing-body interference effects and flow-field calculations have not been 
considered to date. In analyzing the present configuration, a number of 
semi empirical approximations have been employed in estimating the more 
important effects. Calculations for various tail configurations are now 
in progress and these have been utilized to some extent in checking the 
reasonableness of the approximations for the tail-body contribution. 


Detailed Considerations and Results 


In the absence of information concerning the end-plate effects of 
bodies and horizontal tails on yawing vertical tails , two methods were 
used to evaluate the vertical-tail contribution to the stability deriva- 
tives Cy r , C nr , and Cj . The first method was an approximation based 

on the results obtained for the sideslipping motion and utilized the 
following relationships 



ap a 

2V 



* a 2* 


2V 


(15) 


where the minus sign is introduced to maintain the conventional system of 
positive forces and moments and 

p _ rl _ rb l 
V 2 V b/2 

It is clear then that 

dg _ _l_ 

9 £b b/2 

2V 


and hence 



(16) 
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Analagous procedures yield 


and 


Wt = "(S)t b/2 

W t “ -(%\ ^ 


( 17 ) 


( 18 ) 


In the preceding equations 2 is the distance between the center 
of gravity of the airplane and the center of pressure of the vertical 
tail (in the presence of the body) for the sideslip motion. 

A second method of estimating the tail contribution to the yawing 
derivatives is to calculate the two limiting cases (similar to what was 
done for sideslip) corresponding to zero- and complete-end-plate situations 
and then obtain an average value. In order to obtain the complete end -plate 
solutions for Cy r and C nr , the results for C L ^ and C m ^ available 

for wings may be utilized by simply making appropriate changes in the 
nondimens ionali zing parameters. The stability derivatives C L ^ and C^ 

are available in reference 1 for subsonic leading edges and in reference 5 
for supersonic leading edges. At present there is no information for the 
zero-end-plate solution corresponding to the subsonic-leading-edge condi- 
tion; for supersonic leading edges the following expressions have been 
obtained: 



5 

KBC 

5A yt 2 (l - K)(KBC - l) 2 (B 2 C 2 - 1 )'|k(BC + 1)\|k(BC - l) 
4KB 5 C 5 (^ - 9K + K^) + 4B 2 C 2 (-2 + 6 k + K 2 - 5K3) + 
ABC(-l - 5K + 10K 2 ) + k(3 - 5K) 



8K 2 B^C^ ( -l + 2K) + 


( 19 ) 
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B c. 


KBC 


Lr ;t,o 


15A 


^3(1 - K) 2 (KBC - 1)3(b2c2 - 1)1k(BC + l)ix(BC - l) 


i4ok 5 b 5 c5(-i + 5K - 5K 2 ) + 


2K 2 B^C^ (6l - 193K + 213K 2 - 21K3) + KB3c 5 (-112 + - 
472K 2 - 3Yk3 + 1381 ^) + B 2 C 2 (30 - 152K + 89 K 2 + 390 K 3 - 
kk9¥^ + 12 K 5 ) + BC(15 + 122K - 1+4 IK 2 + 436k 5 - 12K^) + 

(-43 + 130K - 125K 2 7| (2°) 


Equations (19) and (20) are applicable to swept! ack vertical tails of 
arbitrary aspect ratio and zero taper ratio with supersonic leading 
and trailing edges; rotation and moments are measured about the apex 
of the vertical tail. 


The differences between the two methods of estimating the tail con- 
tribution to the yawing derivatives may be evaluated. Figure 21 presents 


the results for 




t 


based on equation (l6) and also the limiting 


cases corresponding to zero- and complete -end -plate situations using 
equation (19) and wing results for C^. numerical results pre- 

sented in figure 21 are for a yawing motion about the center of gravity; 
hence the usual transfer -of -axis formula, which involves the previously 
mentioned , must be used in conjunction with equation (19).^ 

It is seen that an average value for (pY r )t ^ ase( ^- on th e second method 


would not be significantly different from the value calculated by equa- 
tion (16). A similar situation was found to exist for the derivative 

In view of this evidence, it was not felt necessary to estimate the zero- 
and complete -end -plate solutions for ( c Z r )t* Tilus equation (l8) was used 

to calculate the tail contribution to C i r . 


With regard to the body contribution to the side force and yawing 
moment, the difficulties previously encountered for sideslip motion are 
also present for the yawing motion. These have been treated in an 
analagous manner (see sideslip section) using equivalent bodies of revo- 
lution and results obtainable in references 15 and 1 6. The body contri- 
bution to the rolling moment was assumed to be small and hence was 
neglected. 

The derivatives CY r > C nr , and for the complete configuration 

(and for the body and tail contributions) are presented in figures 22, 

23, and 24, respectively. 
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CONCLUDING REMARKS 


Theoretical calculations of the stability derivatives at supersonic 
speeds for a high-speed airplane configuration have been presented. The 
methods used represent detailed consideration of the important effects 
and contributions of the various airplane components utilizing available 
theories and theoretically justifiable and reasonable approximations. 

The results are valid for motions involving small angles of attack and 
sideslip and low rates of rolling, pitching, and yawing. 

It should be emphasized that the theoretical information available, 
although supplying many of the basic requirements, is by no means suffi- 
cient to allow estimation of the various forces, moments, and stability 
derivatives for arbitrary complete configurations. Some of the approxi- 
mations developed in the present analysis would be applicable; others 
would require considerable modifications depending on the configuration 
under study. It is felt, however, that a detailed analysis analagous to 
that undertaken for the present airplane would result in reliable esti- 
mates for many types of complete airplane configurations suitable for 
supersonic flight. Further theoretical research on fundamental problems 
indicated in the text is of course required in order to improve the 
accuracy of the estimations and to help minimize the time and effort 
required to obtain them. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics , 
Langley Field, Va. , July 27, 1955* 
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TABLE I.- CHARACTERISTICS OF AIRPLANE CONFIGURATION 
[Dimensions correspond to those of a l/lO-scale model] 


Over -all length, in 

Wing area (total), sq. ft 

Wing airfoil section . 

Horizontal-tail area (total), sq. ft . . . 
Horizontal-tail airfoil section . . . . . 
Vertical-tail area (exposed), sq. ft . . . 

Vertical-tail airfoil section 

Center-of-gravity location, percent M.A.C 


.... 90 
. . . 4.01 

NACA 65AOO3 
. . 0.908 
NACA 65AOO3 
. . 0.874 

NACA 65 A003 
.... 52 
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Figure 1. - Continued. 
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Figure 8. - Comparison at Mach number { 2 between the span load distributions 
for the rolling wing-body combination and the corresponding isolated wing 
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Figure 9*- Finite vortex used to approximate a tent lifting line 
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Figure 11.- Induced dovnwash distributions on the horizontal tail for 
rolling motion. Note that negative values for w denote downwash 
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Figure 13*- Procedure used to simulate the induced sidevash (and downwash) 
distributions by means of superposing elementary solutions. 
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Figure 21.- Comparison of the tail contribution to the stability deriva- 
tive CY r calculated from equation (l6) with the complete- and no- 

end-plate solutions. 
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